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Abstract 

Using a confining, Dyson- Schwinger equation model of QCD at finite tem- 
perature (T) and chemical potential ()«) we study the behaviour of (qq), m^, 
/,r and the masses of the p- and (/>-mesons. For each of these quantities there 
is a necessary anticorrelation between its response to T and its response to [i. 
(—(qq)) and f n decrease with T and increase with /z; m v is almost insensitive 
to T and [i until very near the border of the confinement domain; and the 
mass of the longitudinal component of the vector mesons increases with T 
and decreases with //. At T = 0, the p-meson mass is reduced by approxi- 
mately 15% at nuclear matter density. These results are a consequence of the 
necessary momentum-dependence of the dressed-quark self-energy. 
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High-energy heavy-ion (A- A) experiments, which produce systems with large baryon 
density, are an important preliminary step in the search for a quark-gluon plasma. An 
outcome of these experiments is the observation that the dilepton yield in the region below 
the p-resonance is approximately five-times greater than that seen in proton induced (p-A) 
reactions [IJ. One model calculation |2| shows that this enhancement can be explained by a 
medium-induced reduction of the p-meson's mass, another || that it follows from an increase 
in the p-meson's width. A decrease in the p-meson's mass is consistent with the QCD sum 
rules analysis of Ref . [|j] but inconsistent with that of Ref . , which employs a more complex 
phenomenological model for the in-medium spectral density used in matching the two sides 
of the sum rule. In Ref. |J there is no shift in the p-meson mass but a significant increase 
in its width. The consistency between Refs. || and f| is not surprising since, in contrast 
to Ref. [f|] , they both rely heavily on effective Lagrangians with elementary hadron degrees- 
of-freedom. Herein, using a simple, Dyson- Schwinger equation (DSE) model of QCD, we 
reconsider the response of the p-meson's mass to increasing baryon density. In focusing on 
dressed-quark and -gluon degrees of freedom, our study has similarities to that of Ref. j| . 

DSEs provide a nonperturbative, continuum framework for analysing quantum field the- 
ories. While they have been used extensively at T = = p in the study of dynamical 
chiral symmetry breaking (DCSB) and confinement H], and the calculation of hadronic 
observables [[!]], their application at finite T and p, though straightforward, is in its infancy. 

We begin with the specification of a model dressed-gluon propagator 

g 2 D,M n*) = (V - 2vr 3 1 4o 5\v) , (i) 

where (p M ) := (p,Q k ), ^fc = 2/otT is the boson Matsubara frequency and rj is a mass-scale 
parameter. Equation ([!]) is an extension to finite-T of the model for introduced in 
Ref. ||. It has the feature that the infrared enhancement of the dressed-gluon propagator 
suggested by the DSE studies of Refs. || is manifest. However, it represents poorly the 
behaviour of D^ u (p,Q k ) away from |p| 2 + 11^0. Nevertheless, this limitation leads only 
to easily identifiable artefacts and hence does not preclude the judicious use of the model. 
Another feature of (|1|) is that it specifies a model with the merit of simplicity: it provides 
for the reduction of integral DSEs (e.g., the gap equation and Bethe-Salpeter equations) to 
algebraic equations, which facilitates the elucidation of many of the qualitative features of 
more sophisticated models. 

Using (Jlj) and the rainbow-truncation for the dressed-quark-gluon vertex: T a {k,p) = 
|A a 7^, the DSE for the dressed-quark propagator in our model for QCD at finite-T and p 
(QCD£) is P 



S 1 (p, u k ) = S \p, uj k ) + \r} 2 7„S{pi u k )>y u , (2) 

where Sq 1 (p, oj k ) '■= ij -p + ilA UJ k+ + m , with u k+ = u k + zp, u k = (2k + 1)ttT is the fermion 
Matsubara frequency and p is the chemical potential, and m is the current-quark mass. 
(In our Euclidean formulation, {7^,7^} = 25^ u with 7^ = 7^.) The solution of (g) has the 
general form 

k z*7 ■pA(p k ) + i^iU k+ C(p k ) + B(p k ) ' 

= -ij ■ p a A (p fc ) - i^uj k+ a c (p k ) + &b (Pk) , (4) 
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where p k : = (p,u k+ ), and (§) entails that the scalar functions introduced in (||) satisfy 

rfm 2 = B 4 + mB 3 + (Ap 2 k - if - m 2 ) B 2 - m (2 rj 2 + m 2 + 4pfj B , (5) 

A(p k ) = C{p k ) = W{ ^1 . (6) 

In the chiral limit (m = 0) the character of the solution of (§) is transparent. There are 
two qualitatively distinct solutions. The "Nambu-Goldstone" solution, for which 



otherwise 

Re(p 2 k ) < £ 

C ^ = ^ ^l + ,/TT¥^, otherwise, (8) 
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describes a phase of this model in which: 1) chiral symmetry is dynamically broken, because 
one has a nonzero quark mass-function, B(p k ), in the absence of a current-quark mass; and 
2) the dressed-quarks are confined, because the propagator prescribed by these functions 
does not have a Lehmann representation. The alternative "Wigner" solution, for which 



B (p k ) = 0, C (p k ) = I 1 + Jl + |r , (9) 



2 V \ Pfe 

describes a phase of the model in which chiral symmetry is not broken and the dressed-quarks 
are not confined. 

With these two phases, characterised by qualitatively different, momentum-dependent 
modifications of the quark propagator, this model of QCDj. can be used to explore simulta- 
neously both chiral symmetry restoration and deconfmement. That and the calculation of 



its equilibrium thermodynamic properties are the subject of Ref. [II|. The model exhibits 
coincident, first order deconfmement and chiral symmetry restoration for all /x 7^ but the 
coincident transitions are second order for fj, = 0. The extreme points on the phase boundary 
are (T = 0,/i w 0.28 77) and (T w 0.16 77, fi = 0). 

The vacuum quark condensate is proportional to the matrix trace of the chiral-limit 
dressed-quark propagator. Using (^[), (0) and (|j) we obtain the following expression, valid 
in the domain of confinement and DCSB 

- m =V 3 -fTY, dy y 2 Re U\ - y* - u 2 + ) , (10) 

7T ;=Q JO \ / 

with: f = T/r], = fi/ri, u) i+ = uJi+/v, and uf maK <\+fi 2 , Af = \ + Ji 2 -u) 2 . For T = = fi, 
(— (qq)) — ?7 3 / (80 7r 2 ) = (0.11 ?7) 3 . In Fig. [I] we see that (— (qq)) decreases with T but 
increases with increasing /1, up to a critical value of ix c {T) when it drops discontinuously to 
zero.Q The increase of (—(qq)) with /1 must be expected in the confinement domain because 



1 These results are in qualitative and semiquantitative agreement with (T = 0, fi 7^ 0) |T^] and 
(T / 0,/i = 0) |L3| studies of a more sophisticated model that better represents the behaviour 
of in the ultraviolet. The /^-dependence is also qualitatively identical to that observed in a 
random matrix theory with the global symmetries of the QCD partition function [14|. 
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confinement entails that each additional quark must be locally-paired with an antiquark, 
thereby increasing the density of condensate pairs. This vacuum rearrangement is manifest 
in the behaviour of the necessarily-momentum-dependent scalar part of the quark self energy, 
B(p k ). 

Our primary interests are the bound state properties of tt- and p-mesons. The tt has 
been much studied and it follows [[HJ from the axial- vector Ward-Takahashi identity that, 
writing the it Bethe-Salpeter amplitude as 



T n (p; P)=i l5 B (p 2 



(11) 



one can obtain a reliable approximation in the calculation of those it observables for which 
the dominant, intrinsic momentum- scale is less-than 10 GeV 2 . Using (|ll|) then 0)0] i n the 
domain of confinement and DCSB 



^ j y "max P^H s / 

(m qq) w = rf — ^ / dyy 2 Re [B [a Bo - B 

TT ;=Q JO 



-2 -2 , 2-2 , -2 

u l+ a c + y a A + a B 



)} 



(12) 
(13) 



where \p\ — rjy and Bo(fti) := r)Bo(j)y, TjOi), etc. The right-hand-side of ( pT3| ) is zero for 
m = and increases linearly with m, for small-m. In (^) the canonical normalisation 
constant for the it Bethe-Salpeter amplitude is 



O AT t max 



TT 



1=0 



dy y 2 Re [B 2 [a 2 A - 2 [of + a c a' c + y 2 o A a' A + a B a' B ] (14) 



-h 2 



UJ 



i+ 



ac^-(a' c ) 2 )+y 2 (a A a'i-(a 



r' A f 



-J ^2 



with a' := dcr(y,ui + )/dy, which provides a quantitatively accurate approximation to the 
leptonic decay constant .0 

In the chiral limit, we have from that 



B 



BV 



(15) 



oc = C = a a and a' c = = a' A . The simplicity of the model is again manifest in these 
identities, which yield 



fl = V 



lQN r - 



In 



71- 



Af 



l + 4/i 2 -4^-|An . 



-,2 8 A 2 



(16) 



z=o 



Characteristic in ( |i~6|) is the combination \j? — uf, which entails that, whatever change f n 
undergoes as T is increased, the opposite occurs as /i is increased. Without calculation, 



2 The relation between the normalisation of the tt Bethe-Salpeter amplitude and the leptonic decay 
constant is discussed in Ref. [15|. The demonstration that (|i~4|) provides an accurate estimate 



of the tt decay constant when using (11) is an antecedent to Ref. [[LJ 
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(|16|) indicates that f n will decrease with T and increase with fi, and this provides a simple 



elucidation of the calculated results in Refs. [|12] , P~3|| . Figure |2] illustrates this behaviour for 

In Fig. H] we also plot m n , from ([12]). It is insensitive to changes in /z and only increases 
slowly with T, until T is very near the critical temperature. This insensitivity is the result of 
mutually cancelling increases in (mqq) v and f^, and is a feature of studies that preserve the 
momentum-dependence of the confined, dressed-quark degrees of freedom in meson bound 
states. 

With 7] = 1.37 GeV and m = 30MeV, one obtains f n = 92MeV and = 140 MeV 
at T — = fi. That large values of rj and m are required is a quantitative consequence 
of the inadequacy of ([!]) in the ultraviolet: the large-p 2 behaviour of the scalar part of the 
dressed-quark self-energy is incorrect. This defect is remedied easily [O] without qualitative 



changes to the results presented here |17 



With the vector Ward identity unable to assist with a significant simplification of the 
bound state problem, p-meson properties are more difficult to study: one must solve di- 
rectly the vector-meson Bethe-Salpeter equation. The ladder truncation of the kernel in the 
inhomogeneous axial-vector vertex equation and the rainbow truncation of the quark DSE 
form an axial- vector Ward-Takahashi identity preserving pair [|10|. It follows that the ladder 



BSE is accurate for flavour-nonsinglet pseudoscalar and vector bound states of equal-mass 
quarks because of a cancellation in these channels between diagrams of higher order in the 
skeleton expansion of which this pair of truncations is the lowest order term. 

A ladder BSE using the T = limit of (|l|) was introduced in Ref. ||. It has one notable 
pathology: the bound state mass is determined only upon the additional specification that 
the constituents have zero relative momentum. This specification leads to a conflict with 
the axial- vector Ward-Takahashi identity, which relates the momentum dependence of the n 
Bethe-Salpeter amplitude to the functions in the quark propagator |I5| . We find this to be 
an artefact of implementing the delta-function limit discontinuously; i.e., the identities [T5| 
between the scalar functions in the 7r Bethe-Salpeter amplitude, and A(p 2 ) and B(p 2 ) are 
manifest for any finite-width representation of the delta-function, as this width is reduced 
continuously to zero. In other respects this ladder BSE provides a useful qualitative and 
semi-quantitative tool for analysing features of the pseudoscalar and vector meson masses. 
For example, Goldstone's theorem is manifest, in that the 7r is massless in the chiral limit, 
and also m\ rises linearly with the current-quark mass. Further, there is a naturally large 
splitting between m n and m p , which decreases slowly with the current-quark mass. 

To illustrate this and determine the response of m p to increasing T and /i, we generalise 
the BSE of Ref. § to finite- (T, fi) as 

T M (p k ; Pe) = -j Re { 7 „ S(p { + \P t ) T M {p t] P t ) Sfe - \P t ) 7/i } , (17) 

where Pi := (P,Qg). The bound state mass is obtained by considering Pi =0 and, in ladder 
truncation, the p- and cj-mesons are degenerate. 

As a test we first consider the n equation, which admits the solution 

^(Po) = 75(^1 + 7-^2) (is) 
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and yields the mass plotted in Fig. |3|. The mass behaves in qualitatively the same manner 
as m n in Fig. ||, from (12|), as required if ( |17|) is to provide a reliable guide. In particular, it 



vanishes in the chiral limit. 

In the case of the p-meson there are two components: longitudinal and transverse to P. 
The BSE has a solution of the form 



r p= {(*_ i Bz.An > (19) 




where 6 P+ labels the longitudinal and 6 P _ the transverse solution. The eigenvalue equation 
obtained from (|T7|) for the bound state mass, M p ±, is 

£ Re {<x 5 « - ±M 2 ± ) 2 - [± < - lM p 2 ± ] a y « - i^ ± ) 2 } = 1 . (20) 

The equation for the transverse component is obtained with [— Uq + — \M^_] in (PUD- 
Using the chiral-limit identities, (|l"5|), one obtains immediately that 



M 2 _ = i?] 2 , independent of T and /i. (21) 

This is the T = = \x result of Ref. ||. Even for nonzero current-quark mass, M p _ 
changes by less than 1% as T and \l are increased from zero toward their critical values. 
Its insensitivity is consistent with the absence of a constant mass-shift in the transverse 
polarisation tensor for a gauge-boson. 

For the longitudinal component one obtains in the chiral limit: 

M 2 p+ = \if - 4(/i 2 - tt 2 T 2 ) . (22) 

The characteristic combination [pi 2 — 7r 2 T 2 ] again indicates the anticorrelation between the 
response of M p+ to T and its response to /i, and, like a gauge-boson Debye mass, that M 2 + 
rises linearly with T 2 for fi = 0. The m^fl solution of ( ^0|) for the longitudinal component 
is plotted in Fig. [| As signalled by (p2|), M p+ increases with increasing T and decreases as 
H increases.^ 

We have stated that contributions from skeleton diagrams not included in the ladder 
truncation of the vector meson BSE do not alter the calculated mass significantly because 
of cancellations between these higher order terms [ffD]. This is illustrated explicitly in two 



calculations: Ref. ||18|| , which shows that the p — > txtx — > p contribution to the real part of 
the p self-energy; i.e., the tt-tt induced mass-shift, is only —3%; and Ref. jn|, which shows, 
for example, that the contribution to the cu-meson mass of the u — > 37r-loop is negligible. 
Therefore, ignoring such contributions does not introduce uncertainty into estimates of the 
vector meson mass based on (IT 



3 There is a 25% difference between the value of n required to obtain the T = = p values of m n 
and f n , from (|l3|) and (|l4|), and that required to give M p ± = 0.77 GeV. This is a measure of the 
quantitative accuracy of our algebraic model. 
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Equation ( |20| ) can also be applied to the 0-meson. The transverse component is insensi- 
tive to T and /i, and the behaviour of the longitudinal mass, M^ + , is qualitatively the same 
as that of the p-meson: it increases with T and decreases with /i. Using 77 = 1.06 GeV, the 
model yields M 0± = 1.02 GeV for m s = 180 MeV at T = = fi. 

In a 2-flavour, free-quark gas at T = the baryon number density is ps = 2/i 3 /(37r 2 ) , by 
which gauge nuclear matter density, po = 0.16 fm~ 3 , corresponds to fi = /io := 260 MeV = 
0.245 7]. At this chemical potential our model yields 

M p+ (/i ) w 0.75M p+ (/i = 0) , (23) 
M^+(/i ) ~ 0.85M^+(/i = 0) . (24) 



The study of Ref. |L2| indicates that a better representation of the ultraviolet behaviour 
of the dressed-gluon propagator expands the horizontal scale in Fig. with the critical 
chemical potential increased by 25%. Based on this we judge that a more realistic estimate 
is obtained by evaluating the mass at fi' = 0.20 rj, which yields 

M P M) «0.85M p+ (/i = 0), (25) 
« 0.90M^ + (/i = 0) ; (26) 



a small, quantitative modification. The difference between (E3) and fl25|), and that between 



(23j) and (P6|), is a measure of the theoretical uncertainty in our estimates in each case. 
This reduction in the vector meson masses is quantitatively consistent with that calculated 
in Ref. M and conjectured in Ref. PD"| . At the critical chemical potential for T = 0, 
M p+ « 0.65M p+ (p = 0) and « 0.80 M^/i = 0). 

We have analysed a simple, DSE model of QCDj. that preserves the momentum- 
dependence of gluon and quark dressing, which is an important qualitative feature of more 
sophisticated studies. The simplicity of the model means that many of its consequences can 
be demonstrated algebraically. For example, it elucidates the origin of an anticorrelation, 
found for a range of quantities, between their response to increasing T and that to increasing 



fx, discovered in contrasting the studies in Refs. |12| and [13| . 

We find that both (—(qq)) and f n decrease with T and increase with p, and this ensures 
that m v is insensitive to increasing fi and/or T until very near the edge of the domain 
of confinement and DCSB. The mass of the transverse component of the vector meson is 
insensitive to T and fi while the mass of the longitudinal component increases with increasing 
T but decreases with increasing \i. This behaviour is opposite to that observed for (—(qq)) 
and f w , and hence the scaling law conjectured in Ref. |21J] is inconsistent with our calculation, 
as it is with others of this type. 

Our study has two primary limitations. First, we cannot calculate the width of the 
vector mesons in this model because the solution of ( |T7j ) does not provide a realistic Bethe- 
Salpeter amplitude. We are currently working to overcome this limitation. Second, the 
reliable calculation of meson-photon observables at T = = \i only became possible with 



the determination 0] of the form of the dressed-quark-photon vertex. The generalisation 
of this vertex to nonzero T and p is a necessary precursor to the study of these processes at 
T^O^p. 
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FIG. 1. The quark condensate, (10), as a function of \i for a range of values of T. In all existing 
studies, in which the quark mass function has a realistic momentum dependence, it increases with 
(j, and decreases with T. At the critical chemical potential, /x c (T), (—(qq)) drops discontinuously 
to zero, as expected of a first-order transition. For fj, = it falls continuously to zero, exhibiting a 
second-order transition at T c (/jl = 0) = 0.1677. 
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FIG. 2. The pion mass, (|l2"|), and weak decay constant, (|T4|), as a function of fj, for a range 
of values of T. m n falls slowly and uniformly with \x [m n (T = 0,^ c ) = 0.95 m n (T = 0, = 0)] 
but increases with T. Such a decrease is imperceptible if the ordinate has the range in Fig. |3|. f n 
increases with \i and decreases with T [f n (T = 0, /i c ) = 1.51 f n {T = 0, fj, = 0)]. 
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FIG. 3. M p+ and as a function of fi for T = 0,0.1. On the scale of this figure, m n is 
insensitive to this variation of T. The current-quark mass is to = 0.011 rj, which for r] = 1.06 GeV 
yields M p+ = 770 MeV and m w = 140 MeV at T = = fi. 
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